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& • Abstract 

Wc present the derivation of the hydrodynamic limit under Eule- 
rian scaling for a general class of one-dimensional interacting particle 
s ! ' systems with two or more conservation laws. Following Yau's relative 

entropy method it turns out that in case of more than one conservation 
laws, in order that the system exhibit hydrodynamic behaviour, some 
^-j^ . particular identities reminiscent of Onsager's reciprocity relations must 

f^S ' hold. We check validity of these identities for a wide class of models. 

It also follows that, as a general rule, the equilibrium thermodynamic 
entropy (as function of the densities of the conserved variables) is a 
globally convex Lax entropy of the hyperbolic systems of conservation 
laws arising as hydrodynamic limit. The Onsager relations arising in 
. this context and its consequences seem to be novel. As concrete ex- 

£h amples we also present a number of models modeling deposition (or 

domain growth) phenomena. 

> 
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■ 1 Introduction 

a, 

We investigate the hydrodynamic behaviour of a very general class of one 
dimensional interacting particle systems with two or more conserved ob- 
servables. The systems are not reversible and the hydrodynamic limit under 
Eulerian scaling is investigated. We apply Yau's relative entropy method 
and obtain validity of the hydrodynamic limit up to the occurence of the 
first shock wave in the solution of the limiting pde. There is no novelty in 
the standard steps of the relative entropy proof, so we only sketch these. 
The real novelty appears when it turns out that, in case of more than one 
conserved quantity, in order to complete the relative entropy proof, a class 
of identities should hold, relating the macroscopic fluxes appearing in the 
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hydrodynamic pdes. These identities are much reminiscent of Onsager's reci- 
procity relations. As far as we know these relations have not been pointed 
out in the context mathematically rigorous Eulerian hydrodynamics. We 
check the validity of these relations for the very general class of interacting 
particle systems considered. As a consequence of the Onsager relations it 
follows that the systems of partial differential equations (systems of conser- 
vation laws) arising as hydrodynamic limit are by force of hyperbolic type 
and the equilibrium thermodynamic entropy of the system (as function of 
the densities of the conserved quantities) is globally convex Lax entropy of 
the hydrodynamic equations. This fact may be not so surprising, as it is 
commonly accepted physical fact. So much so that hyperbolic systems of 
conservation laws possessing a globally convex Lax entropy are commonly 
called of physical type, see §. Nevertheless, as far as we know, this fact has 
not been emphasized in the context of mathematically rigorous derivation 
of hydrodynamic behaviour. It is worth noting that given a hyperbolic sys- 
tem of conservation laws the existence of convex Lax entropies is far from 
trivial: in the case of two component systems the local existence of con- 
vex Lax entropies was established in the very technical work 0. In case of 
more than two components in general the pdes defining Lax entropies are 
overdetermined, so in general Lax entropies do not exist at all. It turns out 
from our result that the hyperbolic systems of conservation laws arising as 
hydrodynamic limit are of very special type: they always possess a globally 
convex Lax entropy, namely the equilibrium thermodynamic entropy of the 
system. 

Beside the general framework we also present a number of concrete ex- 
amples of deposition models with two conserved quantities to which the 
general result applies, deriving in this way systems of pdes (hyperbolic sys- 
tems of conservation laws) which describe macroscopically domain growth 
phenomena in 1+1 dimension. 

Our general results are easily extended to more than one dimensions, 
however the formalism becomes more complicated. We were mostly mo- 
tivated by the (one dimensional) deposition models presented as concrete 
examples. 

The paper is organized as follows: In section ^ we present the general 
formalism and the conditions under which the hydrodynamic limit is derived. 
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In section ^ we state the main results of the paper. In section || we present 
a number of concrete examples to which the general framework applies. We 
hope that the models introduced in this section could be of interest in the 
context of deposition/domain growth phenomena. In section || we sketch 
the proof of the main result formulated in section ||| The sketchy proof is 
broken up into several parts. We only hint at the standard steps of the 
relative entropy proof, referring the reader to the original work [^] or the 
monographs § or |{. The essential parts of this section are subsections 



.2| and 1^31 where the Onsager relations and their consequences are derived. 
Finally, in section || we extend the results formulated in the previous sections 
from two to arbitrary number of conserved quantities. 



2 Microscopic models 

2.1 State space, conserved quantities 

Throughout this paper we denote by the discrete tori Z/iVZ, N € N, 
and by T the continuous torus M/Z. We will denote the local spin state 
by S, we only consider the case when S is finite. The state space of the 
interacting particle system is 

Q N := s jN . 

Configurations will be denoted 

For sake of simplicity we consider discrete (integer valued) conserved quan- 
tities only. The two conserved quantities are denoted by 

£ : S -> Z, 
T) : S -> Z, 

we also use the notations £j = £(ct>j), rjj = r](u)j). This means that the 
sums £j and ^ ■ rjj are conserved by the dynamics. We assume that the 
conserved quantities are different and non-trivial, i.e. the functions £,77 and 
the constant function 1 on S are linearly independent. 
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2.2 Rate function, infinitesimal generator 

We consider the rate function r : SxSxSxS^ R+. The dynamics of 
the system consists of elementary jumps effecting nearest neighbour spins, 
(uj,ujj+i) — ► (<jj'pUj +1 ), performed with rate r(u)j,Uj+i;Uj } u)j +1 ). 
We require that the rate function r satisfy the following conditions. 



(A) If r(ui,U2;co' 1 ,Lj' 2 ) > then 

7/(wi) +ri(u 2 ) = n^) +n(uJ 2 ). 

(B) For every K G [min£, max£], L G [min 77, max 77] the set 



(1) 



K,L 



j£T N j£T N 



is an irreducible component of O , i.e. if u,u' G then there exists 
a series of elementary jumps with positive rates transforming u into a/. 

(C) There exists a probability measure tt on S such that for any u±, cj 2 , uj[, 

7r(cui)7r(cJ2)r(cJi,cJ 2 ;cj / 1 ,u;2) = 7r(u;07r(c4)K w 2> w i; w 2,^l)- 

(D) We denote R(uji,uj2) ■= Y^u}' 1 ,u)' 2 eS r ( u "L> UJ 2> UJ ii UJ 2)- This is the total 
jump rate of the nearest neighbour spin pair (wi,w 2 )- Then for any 

R(ui,U)2) + i?(w2,W 3 ) + ii(w 3 , Wi) = i?(wi,CJ 3 ) + R(t0 3 ,UJ2) + R(u>2,Vl)- 

For a precise formulation of the infinitesimal generator on Q N we first 
define the map 6j V " : ft N -> 0* for every G S, j G T^: 

cj' if i = j 
cj" if i = j 

if + 

The infinitesimal generator of the process defined on Q N is 

L N f(u)=J2 E K^^i + i;^^")(/(©" v "^)-/(^))- 

jgT^ u',u"eS 



&j u L )={u)" if i= j + 1 
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We denote by the Markov process on the state space Q N with infinites- 
imal generator L N . 

Remarks: 

(1) Condition ([A]) implies that Ylj Cj an d Ylj Vj are indeed conserved quan- 
tities of the dynamics, while condition (|b]) ensures that there are no 
other hidden conservation laws. 

(2) Condition (^) is somewhat imlicit. It seems to us that it is far not 
trivial (if not impossible) to formulate explicit conditions involving the 
rate functions which would be necessary and sufficient for irreducibility. 
However, in the concrete examples treated in section |] one can easily 
check that irreducibility holds. 

(3) Conditions ([A]), (||) and (|C]) determine the measure 7t(lo) up to an ex- 
ponential distortion, that is the probability measures satisfying these 
conditions are of the form @ of the next subsection. 

(4) Conditions (|C|) and (|D|) imply that the stationary measures of the pro- 
cess Xf are computable and have the structure required for hydrody- 
namic behaviour. See the next subsection for details. Another conse- 
quence of these conditions is Lemma |l] which turns out to be of crucial 
importance for hydrodynamic behaviour. 

2.3 Stationary measures, reversed process 

For every S,t E R let G(9,t) be the moment generating function defined 
below: 

G(9, t) := log e e ^ +Tr >^ir(uj). 

In thermodynamic terms G(9,t) corresponds to the Gibbs free energy, see 
JjJ. We define the probability measures 

vr 0iT (u;) := tt(w) exp(0£(w) + 777(0;) - G(9, r)) (2) 

on S. 
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Using conditions (Q) and (|D|), by very similar considerations as in Jl[], 
||] or [Q| one can show that for any 6 and r the product measure 

is stationary for the Markov process on on £l N with infinitesimal gener- 
ator L . We will refer to these measures as the canonical measures. Since 
Y2j £j an d are conserved the canonical measures on Q N are not er- 

godic. The conditioned measures defined on Q, 1 ^ L by: 



N r \ N I 



E& = ^E^ = L 



are also stationary and due to condtion ([B]) satisfied by the rate functions 
they are also ergodic. We shall call these measures the microcanonical mea- 
sures of our system. (It is easy to see that the measure 7r^ l does not depend 
on the values of 6, r.) 

The elementary movements of the reversed stationary process are (wj-i, oJj) 
— ► (w'j-iiWj) with rate r(ujj, cjj-i] u/-, k>j_i)- The reversed generator is 

jgT^ u',ul"eS 

This is the adjoint of the operator L N with respect to all microcanonical 
(and canonical) measures. I.e. the reversed process is the same for any ttq t 

N 

K,L- 



N 

or 7T 



2.4 Expectations 

Expectation, variance, covariance with respect to the measures ttq will be 
denoted by E 9)T (.), Var 0jT (.), Cov e>r (.). 

We compute the expectations of the conserved quantities with respect 
to the canonical measures, as functions of the parameters 9 and r: 

u(6,t) := E, ir (e) = E^KrM = ^(0,r), 

ujes 

v(6,t) := B etT ( V ) = Y,v(^eA^) = G' T (^r). 
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Elementary calculations show, that the matrix-valued function 

U <)-\G» T G^J-Gfrr) 

is equal to the covariance matrix CoV0 )T (£, rj) and therefore it is strictly pos- 
itive definit. It follows that the function (6,t) i— > (u(9,t),v(9,t)) is invert- 
ible. We denote the inverse function by (u,v) t— > (9(u,v),t(u,v)). Actually, 
denoting by (u, v) i— > S(u, v) the convex conjugate (Legendre transform) of 
the strictly convex function (6,t) h- > G(9,r): 

S(u,v) := sup (u6 + vt-G(9,t)), (3) 

we have 

0(u,v) = S' u (u,v), t(u,v) = S' v (u,v). 

In probabilistic terms: S(u, v) is the rate function for joint large deviations 
°f (Ylj€jiYlj r lj)- I n thermodynamic terms: S(u,v) corresponds to the 
equilibrium thermodynamic entropy, see Q. Let 

(5 5) = fS S ™)=:S"(u,v). 

It is obvious that the matrices G"(6,t) and S"(u,v) are strictly positive 
definit and are inverse of each other: 

G"(e,r)S"(u,v)=I, (4) 

where either (0,t) = (u(6,t),v(9,t)) or (u,v) = (9(u,v),t(u,v)). With 
slight abuse of notation we shall denote: -Ke(u,v),r{u,v) ='■ ^u,v, K(u,v),t(u,v) =: 

^u,vi ~^9(u,v),t(u,v) =: E Ujl) , etc. 

We introduce the flux of the conserved quantities. The infinitesimal 
generator L N acts on the conserved quantities as follows: 

L N ii = -<f>(u)i,Wi+i) + <j)(ui-i,Ui) =: -<j>i + 4>i-i, 
L N r]i = -ip(uji,uJi + i) + ip(iVi-i,iVi) =: -tpi + 
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(5) 



where 

0(wi,w 2 ) := ^2 r(u 1} u 2 ;uj[, w 2 )(£(w 2 ) - £(w 2 )) 

= r ( a 'i' a; 2;wi,w 2 )(C(^i) -C(^i)), 

^(^i,cj 2 ) := ^ r(wi,o; 2 ;a;i,a; 2 )(7?(w2) -??(^2)) 

= r ( w i' a; 2;wi,w 2 )(^(wi) -^i))- 

We shall denote the expectations of these functions with respect to the 
canonical measure v by 

*(«,«) := E U)W (0) 

= ^2 r ( w i> w 2;^i,w 2 )(C(^ 2 ) - C(^2))7r u ,u(wi)7ru,u(^2), 

(6) 

:= E Wj „(V>) 



The first derivative matrix of the fluxes and ^ will be denoted 

°<-) < 7 > 

As a general convention, if 5 : S m — > M. is a local function then its 
expectation with respect to the canonical measure 7r™ i; is denoted by 

A(u,v) := E U) „(<5) = ^ <J(wi,...,o; m )7r W)U (a;i)---7r M)1 ,(a; m ). 

o;i,...,£j ro e5 m 

3 The hydro dynamic limit 

We will show, applying Yau's relative entropy method, that under Eulerian 
scaling the local densities of the conserved quantities u(t,x), v(t,x) evolve 
according to the following system of partial differential equations: 

d t u + d x $(u,v) = , , 

d t v + d x ^(u,v) = 0. 1 j 

It also turns out from our proof (more precisely as a consequence of the 
Onsager relations proved in Lemma fil) that the systems of conservation 
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laws (§) arising as hydrodynamic limits are necessarily of hyperbolic type 
and the equilibrium thermodynamic entropy function (u, v) i— > S(u, v) is a 
(very special) globally convex Lax entropy for the system (|8|). (See H or [ll]] 
for the pde notions used.) This may be not so surprising, as it is commonly 
accepted fact and drops out automatically, without any computations in 
some particular model systems investigated so far. Nevertheless, we have 
not found a general statement or proof of this fact in the hydrodynamic 
limit literature. 

3.1 Notations 

For the proper formulation of our results we need some more notations. Let 
u(t,x),v(t,x),t € [0, T],x £ S be a smooth solution of (§) (more precisely: 
let it be twice continuously differentiable in both variables). We shall use 
the notations 

6(t,x) := 6(u(t,x),v(t,x)) 
r(t,x) := T(u(t,x),v(t,x)). 

The true distribution of the Markov process X^ 1 at macroscopic time t, 
i.e., at microscopic time Nt is 

:=rfe*p{NtL»}. (9) 

The true distribution will be compared to the following time dependent 
reference measure (also called local equilibrium) on Q, N : 

v t ■= II n u(t4)At,ir (10) 

This measure is not stationary (unless u and v are constant), and the lo- 
cal densities of the conserved quantities are discrete approximations of the 
functions u(t,x),v(t,x). 

We shall use a stationary measure ir N := ttq on £l N as an absolute 
reference measure. The Radom-Nikodym derivatives of the true distribu- 
tion and the time dependent reference measure with respect to the absolute 
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referencee measure are denoted as follows: 
K ■= ^(uL) = eMNtL* N }h%. 



f" := jfrW (ID 
= J] expftfaWt, j-) + T)("Mt, j-) - G(9(t, j-),r(t, j-))} 



3.2 The main result 



Our aim is to prove that if fj, is close to v$ in the sense of relative entropy, 
then fj,f stays close to v t in the same sense uniformly for t 6 [0,T]. If we 
consider two different pairs of smooth solutions (ui(t,x),Vi(t,x)), i = 1,2 
of (Q) it is a simple exercise to show that the relative entropy of the two 
time dependent reference measures is of order x N. This suggests that one 
should prove 

H»(t):=H(tf\u») = o(N), (12) 
uniformly for t £ [0, T]. 

Theorem. Consider an interacting particles system model defined as in 
the previous section which satisfies conditions (^), ^), ftQ) and ^). Let 
$>(u,v) and *&(u,v) be defined as in (j^). 

(i) The system of conservation laws ^) is hyperbolic in the domain (u,v) S 
(min£, max£) x (minry, max?]). Furthermore, the equilibrium thermody- 
namic entropy (u, v) i— > S(u, v) is a globally convex Lax entropy for the 
system (Qj. 

(ii) Let [0, T] x T 3 (t,x) i— ► (u(t,x),v(t,x)) be a smooth solution of 
and let ^ and be the measures defined in respectively, (1C). Then, 
if 

H(^\tt n ) = 0(N), 
and (T/lj holds for t = then it will hold uniformly for t 6 [0, T] . 



Remark: Part (i) of the Theorem is commonly accepted fact. In some 
particular models investigated it simply drops out without any computation. 
However, we do not know about any explicit formulation (or proof) of the 
general fact stated here. 
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Prom part (ii) of the Theorem, by applying the entropy inequality in a 
standard way (comparing the true measure ^ with the local equilibrium 
reference measure vf) one gets the following corollary: 

Corollary. Under the conditions of the Theorem, for any t £ [0, T\, the 

following limits hold as N — > oo: 

(i) For any smooth test function g : T — > K. 



jr E 9(j/N)C 3 (t)^ [ g(x)u(t,x)dx, 
jr Yl 9(j/N)vj(t)^ / g{x)v{t,x)dx. 



(ii) The asymptotics of the relative entropy of the true distribution fi t with 
respect to the absolute reference measure tt^ 0>Vq is 

N-^Qj? |<,„ ) [ (S(u(t,x),v(t,x))-S(u ,v ))dx, (13) 

where S(u, v) is the thermodynamic entropy defined in ^). 

Remark: Note that since S(u, v) is Lax entropy of the pde (|8]) the right 
hand side of ([13]) does not change in time as long as the solution (u(t, x),v(t, x)) 
of (||) is smooth, and starts to decrease when the first shock appears. This 
means that the relative entropy H(fj,f \ ^u ,v ) decreases by o(N) before the 
appearence of the first shock in the system. 

4 Examples: deposition models 

If we fix the size of the spin space, then we have finitely many equations from 
the conditions on the rate function, thus we can get a finite-parameter family 
of models. The smallest value of \S\, for which there exists a proper model 
is 3, since we need to have two different non-trivial conserved quantities. 
We present two concrete examples: one with |5| = 3, one with \S\ = 4. 
A third example with \S\ = oo, to which the Theorem applies with some 



modifications, is described in pl|. 



Our concrete examples are deposition models, r/ : 5 — ► N, and £ : S — > Z. 
rjj, respectively, £j are interpreted as particle occupation number, respec- 
tively, (negative) discrete gradient of deposition height. The dynamical 
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driving mechanism is such that 

(i) The deposition height growth is influenced by the local particle density. 
Typically: growth is enhanced by higher particle densities. 

(ii) The particle motion is itself influenced by the deposition profile. Typ- 
ically: particles are pushed in the direction of the negative gradient of the 
deposition height. 

It is natural to assume left-right symmetry of the models. This is realized 
in the following way. There is an involution 

R:S^S, RoR = Id 

which acts on the conserved quantities and the jump rates as follows: 

r){Rw) = ri(u), ew = », 

r(Rw2, Rut; Raj 2 , Ru[) = r(u\,U2;uJi, uj' 2 ). (14) 

Correspondingly on the macroscopic level we shall use the traditional 
notation p(t,x) (instead of v(t,x) of the general formulation) and u(t,x). 
The limiting partial differential equations will be also invariant under the 
left-right reflection symmetry: (p(t,x),u(t,x)) i— > (p(t, — x), — u(t, — x)) . 

4.1 A model with \S\ = 3 

The state space is S = {—1,0,1}. The left-right reflection symmetry is 
implemented by R : S — > S, Rlo = -co. The two conserved quantities 
are = u (the spin itself) and t](uj) = 1 — \u\ (the number of zeros). 
(It is easy to see that up to linear combinations these are the only two 
conserved quantities we can define on S.) From condition (|A|) it follows 
that r(u>i,U2',(jJi,(jJ2) > only if u[ = u>2 and uj' 2 = W\. I.e. the dynamics 
comsists of exchanges of nearest neighbour spins. It follows that, without 
any restriction on the rates, condition (^) is satisfied with any probability 
measure tt on S. Our natural parametrization is 

1 — p rb U 

Kp,u{°) = Pi K p , u {±l) = , 

with the parameter range {(p, u) : p G [0, 1], u G [—1, 1], p + \u\ < 1}. 
Condition (Q) is fulfilled if and only if 

r(l,-l;-l,l) -r(-l,l;l,-l) 

= r(l, 0; 0, 1) - r(0, 1; 1, 0) + r(0, -1; -1, 0) - r(-l, 0; 0, -1) 
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holds. The reflection symmetry condition ( |i~4| ) reads 



r(l,0;0,l) =r(0,-l;-l,0), r(0, 1; 1,0) = r(-l,0;0, -1). 

These conditions leave us with 

r(l,-l;-l,l) = a, r(-l, 1; 1,-1) = 2c + a, 

r(0,-l;-l,0) = b, r(-l,0;0,-l) = c + b, 

r(l, 0,0,1) = b, r(0, 1,1,0) = c + b, 

where a, b > and c > max{— 6, — a/2} are free parameters. Without loss 
of generality we may choose c > (otherwise, rename uj := —u). It is 
easy to check that ondition (|B|) is satisfied iff (a + 2c) (b + c) > 0. We are 
not interested in the c = case, since that defines the reversible process 
which would imply diffusive rather than hyperbolic (Eulerian) scaling. By 
fixing an appropriate time scale we choose c = 1. It is easy to compute the 
microscopic fluxes 4>j an d ifij given by formula (||): 

<f)j = - (uj {{u)j - 1)(1 + - 2au)jUj+i + 26(1 + UjUj + i)), 

ipj = b(uj^ +1 - wj) + -(1 - ujj)(l + + 

The macrosscopic fluxes are computed with formula (|6|. Inserted in (^) this 
leads to the hydrodynamic equation: 



d t p + d x {pu) = 
d t u + d x {p + u 2 ) = 



^ _ n ( 15 ) 



This system is known in the pde community as Leroux's equation. The 
system has some very special features: it belongs to the so-called Temple 
class and it was much investigated. For details see |J. Validity of this pde 
in the hydrodynamic limit up to the occurence of shocks follows from our 
general Theorem. 

Remark: It is an easy exercise to see that all models with \S\ = 3 satisfying 
the general conditions ([A], [B|, [C], |D]), without the extra assumption of left- 
right reflection symmetry, are essentially equivalent, in the sense that in the 
hydrodynamic limit they lead to pde-s which can be transformed to (|l5| ) by 
linear combinations of the functions involved. 
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4.2 A finite bricklayer model 

In the following example we give a finite version of the infinite bricklayers 
model introduced in [11]. Let S = {0, 1} x {—1, 1}. The elements of S will 
be denoted u =: (n,z). Left-right reflection symmetry is implemented as 
R : S — > S, R(n, z) = (n, —z). The conserved quantities are = z and 
rj((jj) = n. Condition (^) leaves twenty (possibly) non-zero rates. Due to 
the left-right reflection symmetry conditions eight pairs of rates are equal. 
Using the notation 



r(u 1 ,u 2 ;uj' 1 ,U2) = r f™ 1 



n, n n n, n 



2 . "1 "2 

« 9 ' z'-i ' z!> ' 



in the following table we list the (possibly) non-zero rates, parametrized by 
twelve nonnegative parameters. 



'0 o.o o^ 

+ ' +' 



1. 1 I s 

+ > +' -, 



1.0 l\ 

' +' -) 



1.0 1 



1.10 



1. 1 0\ 

+' +' +) 

1. 1 0\ 

-' +' -) 

1. 1 0\ 

+' -' +; 

'0 1. 1 0\ 

_, +, +, 
ro 1.1 o 



) 



a, 
c, 
(:■ 



■ r, 
: s, 

X, 

: y, 



(0 0. 

1+' -' -' + 

(1 1.1 1 

1+' -' -' + 



1 0. 1 



)=/. r( 
,£)=P, r( 



10.1 



10.0 



1 0.0 



1 0.0 



1 0.0 



1 0.0 



1 0.0 



1) 
i) 
i) 
i) 
i) 
i) 



-b. 

d, 

e, 
/, 
P, 
q, 
r, 

s, 
x, 

- y- 



All the other jump rates are zero. 
It is easy to check that condition 



imposes 



(16) 



and no other restriction. It also follows that the measures ir PjU are of the 
product form 

1 + zu 



7T, 



(jU (n, z) = (np + (1 - n)(l - p)) ■ 



n 



0,1, z = +,- 



(17) 



with the parameters p 6 (0, 1), u 6 ( — 1, +1). 
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Another straightforward computation shows that condition (^) reads 

c+f+p+y = d+e+q+x ^ 
a+f+q+y=b+e+p+x 

So, we we are left with a nine-parameter family of models. Given the for- 
mulas (|B|) we compute the fluxes <f>j and ipj. Using the conditions (16) and 
(18) eventually we get 

2(j)j = (b - a) + (p - q)(rij + n j+1 ) - (b + a)(z j+1 - Zj) 

+ (a + b - e - f - x - y)(rij + n j+1 )(z j+1 - zj) + (a - b)z j+ iZj 
— {a + b + c + d — 2e — 2f — 2x — 2y)njUj + i(zj + i — Zj) 
-(p - q)(nj + n j+1 )z j+1 Zj 
= —(p + q + r + s + x + y) (rij+i - n 



3) 



+(p - g)(nj + n j+1 )(zj + z j+l ) + (y - x)(n j+1 - nj)(z j+1 - zj) 

-2(p - q)n j n j+ i(z j + z j+1 ) 

-(p + q-r-s-x- y){n j+ i - nj)z j+1 Zj 



The macroscopic fluxes are again explicitly computable. From (||) and (|T 
we get 

$(p,«) = {(p-q)p-(a-b)/2)(l-u 2 ), 
V(p,u) = (p - q)p(l - p)u. 

Without loss of generality we may assume p — q > (otherwise rename the 
microscopic variables hj := rij, zj := —Zj). Further on, p = q leads to 
diffusive rather then hyperbolic (Eulerian) scaling of the particle density, so 
we are interested in the p > q cases only. By setting the appropriate time 
scale we can choose p — q = 1 and denote 7 := 2(p-q) • even tualfy we get 
the system of pdes 

r d t p + d x { P {\ - P )u) = 

{ (19) 

1 ^u + a^p-TXi-u 2 )) =0 

In |^] another family of four-state models with two conserved quantities, 
the so-called two channel traffic models are analyzed. These models also 
satisfy conditions (]A|), (|B|), (^) and (|D|). As a consequence our general 
Theorem is applicable to the two channel traffic models, too. 
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About the relation of our four state deposition models (treated in this 
subsection) and the two channel traffic models treated in ||: Due to the 
different symmetry conditions imposed — we impose the left-right reflec- 
tion symmetry described in the first paragraph of this section, while in || 
symmetry between the two traffic channels is imposed — the two families 
of models do not intersect. The one parameter family of partial differential 
equations derived in || essentially differs form our partial differential equa- 
tions (li~9|). (Actually there is no parameter value for which the two pde-s 
are equivalent.) However, the two families of models show many similarities 
and do have common generalizations. 



5 Sketch of proof 



wc 



The present section is divided into four subsections. In subsection 5.1 
present the first steps of the 'relative entropy method' applied. As there is 
no real novelty in this part, we only list the main steps without the compu- 
tational details which are essentially the same as in the original work 13 j of 
Yau or in Chapter 6. of or in |J. 

It turns out that for a general two (or more) component system some 
identity relating the macroscopic fluxes <3? and ^ is essentially needed for 
completing the proof. These relations are reminiscent of Onsager's reci- 
procity relations of nonequilibrium thermodynamics, see e.g. Chapter 10. D 
of |7]] . However an essential difference is worth noting: while the traditional 
Onsager relations are derived under the condition of reversibility of the mi- 
croscopic dynamics, in our case conditions (|A]), (|C|) and (^) are involved 
which do not imply reversability by any means. Seemingly, these relations 
were not explicitly noted so far in the context of mathematically rigorous 
hydrodynamic limits. This omission is probably explained by the fact that 
in the concrete models investigated so far these identities just droped out 
without any computations. 

In subsection 5^ we prove that under the conditions (|A|), (|C|) and 



these identities hold in general. We shall refer to these identities as Onsager 
relations. It also follows from these identities that the systems of partial 
differential equations arising as hydrodynamic limits under Eulerian scaling 
are indeed of hyperbolic type and the thermodynamic equilibrium entropy 
S(u, v) is globally convex Lax entropy of the hydrodynamic equations, as it 
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is commonly assumed. In subsection 5.3 we formulate the consequences of 



the Onsager relations which are crucial for the further steps of the proof of 
the hydrodynamic limit. 



Finally, in subsection 5.4 we sketch the last steps of the proof. Here 
again we follow the standard steps of the relative entropy method, so we 
omit all computational details, referring only to the main stations of the 



proof. For the computational details of subsections 5.1 and 5.4 we refer 
the reader to Chapter 6. of j|] or to [||. However, we warn the reader 
that the omitted details (in particular the last two steps: the control of the 
block replacement and the one-block estimate) are rather sophisticated and 
mathematically deep. 

5.1 First transformations 

In order to obtain the main estimate ( |l~2|) our aim is to get a Grdmwall type 
inequality: we will prove that for every t € [0, T] 

H N (t) - H N (0) < C f H N (s)ds + o(N), (20) 
Jo 

where the error term is uniform in t € [0, T]. Because it is assumed that 
H N (0) = o(N), the Theorem follows. 

For proving ( p0| ) we try to bound (from above) dtH N (s) by const • 
H N (s) + o(N), uniformly for s € [0,T]. We start form the inequality (]2l| ) 
which is derived in [Q] under very general conditions, valid in our case. 

dtH N (t)<N ^Muf-/ ^M^f- (21) 
Jn N It Jn N Jt 

Next we transform the two terms appearing on the right hand side of (|2l|). 



Equation (22) follows form the smoothness of the functions 9(t, x) and r(t, x) 
and from the entropy inequality applied to the measures /if compared with 
the absolute reference measure -n - ^. 

N / —pF-drf = - £ d x 0(t,j/N) / Mtf (22) 
Ja N ft ^ N 

d x r(t,j/N) j fydi4 



N 

H 



+0{l) 
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Equation ( |23[ ) follows from direct computation of the time derivative of the 
function ff . 

-%dtf = *£d t e(t,j/N) (^-u(t,j/N))d^ (23) 
n" ft Jn" 



+ Y, d t r(t,j/N) f Jvi-v(t,j/N))dtf 

■ _ AT J il 



jet 

Next we replace the local variables <pj and tpj in (p2|), respectively, £j and rjj 
in (|23) by their block averages defined as follows: if Sj is a local microscopic 
variable its block average is defined as 

i 5j + --- + 5 j+ i-i 
°i ■ l 

In the following two block-replacements we use again the smoothness of 
the functions 9(t, x) and r(t, x) and the entropy inequality applied to the 
measures /if compared with the absolute reference measure 71"^ . 

- V d x r(t,j/n) / i\)\d[i N t 

j€T N 

+0(1) 



dtft j n 
Jt 



n N Jt 



YdtO(t,j/N)[ (%-u(t,j/N))dtf (25) 
+ Y dtr(t,j/N) f ^-vitJ/Wtf 

■ -nn/V " " 



+o(0 



The last transformation of this first, preparatory part is replacing in (24) the 
block averages <^-, respectively, ^>j- by their equilibrium averages computed 
at the empirical densities, $(£j,m), respectively, ^{Cj^Vj)- The error terms 
appearing in the third and fourth lines of the right hand side of ( |26|) are 



the most important error terms to be controlled by the so called one block 



IS 



estimate towards the end of the proof. 

- *JV fN 



N 



p Tif-JS £JS f 

/ —=P-dtf = - ]T M&3/n) / Hi4Wt (26) 
Jn N ft .^r Jn N 

jeT iv ' /n,r 

- £ ^(t.j/n) I (<t> l 3 -m l A)) dtf 

Before going on with the standard steps of the relative entropy proof we 
need to make a detour. 

5.2 An Onsager type identity 

Lemma 1. Suppose we have a particle system with two conserved quanti- 
ties and rates satisfying conditions §A^), fiQ) and (^). Then there exists a 
potential function (9,t) i— > U(9,t) such that 

$(9,t) :=^(u(9,T),v(9,r))=U' e , 

(27) 

¥(0,t) :=fKfl ) r) I ^,r)) = [/;, 

or, equivalently 

K = %. (28) 



Proof. We prove (Eg). Throughout the forthcoming proof we adopt the 
notation £j := £(u;j), £j := etc - 
From the definitions 

neA^i^A^) = ex P{#(6 + 6) + T(m + m) ~ 2G(6, t)}^^^^), 
and 

n(u 1 )n(u; 2 )e e ^ + ^ +T ^ + ^- 2G( - e ^ {(6 + 6) - 2u(9,t)} = 

V J^Vt J_ £„ _ 9^^ e («i+6+53)+r(m+»)2+r?3) 



^PE^3)K 1+ &-2?3)^ 

(J3GS 
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where Z(6,t) = exp{G(0, r)}. Similarly, 
(7re,r(wi)7r 6 ) ir (a;2)) / T . = 



£ 7r(w 3 )(»7i + r/2 - 2r ?3 )e e ^ + ^+^)+-(^- 
^ ,T) uses 



Hence 



K(0,r) = YWrf S vr(a; 1 )7r( W2 )vr(a;3)e 9 ^ + « 2+ ^) +T ^ +r ' 2+ ^) 
' a;i,a;2,aJ36S 

x (r?i + 772 - 2773) r ( w i^2,Wi,w 2 )(C 2 -6), 

%&t) = J] 7r(a; 1 )7r( W 2)vr(a;3)e 9 ^ + « 2+ « 3 ) +T ^ + ^ + '' 3 ) 

' a;i,a;2,aJ36S 

X + r ( W l' W 2,W 1 ,CJ 2 )(??2 - %), 

For the proof of the lemma it is enough to prove for any K £ [3 min £, 3 max £] 
and L £ [3 min 77, 3 max rj\ 

^ 7r(a)i)7r(o;2)7r(w3)r(a;i ) a;2;wi,W2) 

2 X (^ + % " 2r?3 ^ - 6) - (6 + £2 - 2£ 3 )(77 2 - 772)) = 

From condition ([A|) imposed on the rate functions it follows that in all 
nonzero terms of the above sum one can replace 771+772 by rj' 1 +rj 2 and rj' 2 —r]2 
by 771 —"]][, and similarly for the £-s. After straightforward computations 
this equation (still to be proved) becomes 

Wi,aJ2 1^3,0^ jWjGS: 

t + +m+%L x (A^i,^,^) + A(u, 2 ,c4,"3)) = 

where A : S x S x S ^ Z is defined as follows 

A(u!,uJ2,u 3 ) := £1(772 - ?? 3 ) +£2(773 - 771) +£ 3 (?7l - 

Note that A is antisymmetric regarding permutation of its variables. 

Next, from condition (Q) it follows that the contribution of the two terms 
on the left hand side of the previous equation is the same. Thus, it is to be 
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proved that 

vr(wi)7r(tJ2)vr(c i ;3)r(tJi,tJ2;^,u;2)A(a;i,a;2, w 3 ) = 

uji,W2,lU3,uj' 1 ,uj' 2 £S: 

vi+vt+v'=l ^2 tt(uji)'k(u;2)tt(uj3)R(uj 1 ,U2)A(u 1 ,U2,uJ3) = 

r?i+»?2+»?'=i 

Finally, from the antisymmetry of the function A and condition (|d|) 
imposed on the function R it follows indeed that this last sum equals zero. 

□ 

5.3 Consequences of the Onsager relations 

Relation ( p8| ) is the same as saying that the matrix D(u(9,t),v(9,t)) ■ 
G"(9,t) is symmetric. Using (||) this also reads as 

S"(u,v) ■ D(u,v) = (S"(u,v) ■ D(u,v)y . (29) 

This relation implies that only hyperbolic two-by-two systems of conserva- 
tion laws @ can arise as hydrodynamic limits. Indeed, as the following 
elementary argument shows relation (|2^) can hold with a positive definite 
matrix S" only if D(u,v) can be diagonalized (in the real sense), which is 
exactly the condition of hyperbolicity of the system (||). Indeed, since S" is 
positive definite, we can write 

D = {S")' 1 ' 2 ((S")~ 1/2 (S"D)(S'T 1/2 ) (S") 1 / 2 , (30) 

which means that D is similar to the real symmetric matrix (S") -1 / 2 (S" D)(S") 
and from this the (real) diagonalizability of D follows. Furthermore, (|2^) is 
spelled out as 

S" $' + S" *' = S" <J>' + S" (31) 

which is readily recognized as the the partial differential equation defining 
the Lax entropies of the system (^|). The function F(u, v) := U(9(u, v), t(u, v)) 
is the corresponding (macroscopic) entropy-flux. See || or [[[(]] for the pde 
notions involved. Thus, part (i) of the Theorem is proved. 

Now we turn to two further consequences of Lemma [l] which turn out to 
be of crucial importance in the hydrodynamic behavior. 
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First, the time derivatives of 6 and r are expressed. From the pde (||) it 
follows that 

d t e = -9' u <f>' u d x u - e' u $' v a x v - 9' v v' u d x u - o' v %d x v. 

Using the identity (|3ll) we replace 

0' & = e'& +t'*' -t'*' 

"it V V U ' V U 1 U V 

in the second term of the right hand side. Using also the straightforward 
identities u' T = v' e and 9' v = t' u (see subsection |2.4[ ) we finally get 

d t e = & u a x e + ^' u d x r, (32) 

and by identical considerations 

d t T = & v d x e + m' v d x r. (33) 

Second, due to identity (p7|), 

(d x e(j/N)<S>(u(j/N),v(j/N))+d x T(j/N)y(u(j/N),v(j/N)^ 

j£T N 

= d x U(u(j/N),v(j/N)) = 0(1). (34) 

j£T N 

5.4 End of proof 

Now we return to proving (p0|). Denote 

V§(u, v; u, v) := ®(u, v) - ®(u, v) - &' u (u, v)(u - u) - <&' v (u, v)(v - v) 



and similarly for T>fy(u, v;u,v). Applying (p2|), ( p3| ) and (p4|), from (^) and 
(|26|) we obtain 



ag-^ff ^ = (35) 



n N ft 

d x 6(t,j/n) f N V$(u(t,j/N),v(t,j/Ny,Z l j ,r l l j )dn? 



jer- 

+ £ d x r(t,j/n) f V*(u(t,j/N),v(t,j/Ny,$,rftdii? 

+ ]T 9 x e(t,j/n) I dtf 

+ £ 9 x r(t,j/n) ! fa dtf 



+0(0 
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The first two terms on the right hand side of ([35|) are estimated by the 
entropy inequality, comparing the measure fj, t with the local equilibrium 
measure v t : 



T>$(u(t,j/N),v(tJ/N);$,rfo 
TW(u(tJ/N),v(t,j/N)\Si,Vi) 



(36) 



+ 



d\x t 

< CH N (t) + 0(NI~ 1 ). 



The last two terms in ( p5[ ) are estimated only integrated against time. Ap- 
plying the so-called one block estimate (see e.g. Chapter 5 of H) one gets 



lim lim iV- 1 V ! ds [ 
z^oojv^oo ^? N Jo Jn N 

lim lim iV" 1 V f ds f 

i~,ooN^oo .^fJo Jn N 1 33 



dut = 0, 



d/xf = 0. 



(37) 



This is the only part of the proof where condition (|B]) is used, which ensures 
ergodicity of the Markov process X t on the 'hyperplanes' fi^ L . 

Finally, inserting (^) and p7| ) in ( |35| ) , via we obtain ( ^0[ ) and thus 
the part (ii) of the Theorem is also proved. 



6 Particle systems with several conserved variables 

As noted in the introduction, the results described in the previous sections 
are also valid for particle systems with more than 2 conserved quantities. 

Before we formulate the general results we have to summarize some no- 
tations . Let n > 2 be fixed integer, and £ = (£ , £ 2 , . . . , ) : S — > ffi n the 
vector of conserved quantities. Throughout the present section bold face 
symbols will denote n-vectors. 

We require the rate function to satisfy similar conditions as listed in 
subsection |2^ (in place of conditions ([A]) and (|B|) we need the suitable 
generalizations). For every 6 G M. n we can define momentum generating 
function G(0) as 

G{6) :=log]Te ^Mu;), 

U>£S 
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and the probability measures 



Tr g (u) := 7r(u/)exp(0.£(a/)-G(0)) 

jeT N 

on 5", respectively, on Q N . We define the expectation of the conserved 
quantities with respect to the measure tTq : 

u{6) :=E„(£) = V e G(0). 

One can easily show, that V^G{6) = Covg(u,u) is positive definite. As a 
consequence, the function 6 i— ► u is invertible, and 

e(u)=V tt 5(u), 

where S'(u) is the convex conjugate of G(0): 

S{u) := sup (u • - G(0)). 

We introduce the flux of the vector of the conserved quantities and its ex- 
pectation: 

0(^1, u; 2 ) := ^2 r (. UJ i^ UJ 2;^'i,^2)(^2) ~ C(^2)) 
*(u) := E 0(u) 

= ^ r^i,^;^,^)^^) -4(^2))7ru(^i)7rw(w2)- 

Wl,W 2 , - 

Now we are able to formulate results of the previous sections in the more 
general setting. 

Using the arguments presented in section || one can show that under 
Eulerian scaling the vector of the local densities of the conserved quantities 
u(t, x) evolve according to the following n-component partial differential 
equation: 

d t u + d x &(u) = 0. (38) 

Lemma |l| applies for any two conserved quantities (i ^ j), thus if 

we denote the derivative matrix of the flux vector <&(«) by D(u) := V u $(u) 
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and the second derivative matrix of the thermodynamic entropy S"(u) := 
V 2 S{u), we get 

S"(u) ■ D{u) = (S"(u) ■ D{u)) ] (39) 

Since 5"' (it) is positive definite ( p0| ) implies that D(u) can be diagonal- 
ized which means that the arising system of partial differential equations is 
hyperbolic. Moreover, (|^) spelled out is 

d 2 S $i d 2 S d 2 S $i d 2 S 

H — = H — , (40) 

duidui duj duiduj duj dujdui dui dujduj dui ' 

with 1 < i < j < n. These are exactly the n(n — l)/2 equations defining the 
Lax entropies of the hyperbolic system ( p8| ) . It is well-known that in the case 
of n > 3 only very special ra-component hyperbolic conservation laws possess 
Lax entropies. In general, the defining equations ( |4fj| ) are overdetermined. 
In M these paricular systems of hyperbolic coservation laws are called of 
physical type. From the previous arguments it follows that only physical 
hyperbolic equations can arise as the hydrodynamic limit of an interacting 
particle system satisfying our conditions. 
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